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II. 


ABSTRACT 


NASA IS INTERESTED IN FIELD TESTING FOR POSSIBLE LEAKAGE IN THEIR FUELING 
SYSTEMS. HOWEVER, MANY FUELS ARE HAZARDOUS TO THE EXTENT THAT 


TEST, AND THAT THOSE RESULTS BE CALIBRATED TO SIMULATE THE ACTUAL 
PROCESS. A TECHNIQUE SUCH AS THIS WOULD ALLOW PERSONNEL TO BE ON SITE 
DURING THE TESTING, AND USE TECHNIQUES TO DETERMINE THE BEHAVIOR OF 
THE SYSTEM THAT COULD NOT BE USED OTHERWISE THIS ENDEAVOR ATTEMPTS 


m. 


SUMMARY 


IT WAS DESIRED TO PREDICT THE LEAKAGE OF VARIOUS FLUIDS WITH THE KNOWN 
LEAKAGE OF A KNOWN GAS FOR A FIXED CONFIGURATION. A SIMPLE MODEL WAS 
CONTRIVED, SOFTWARE WAS DEVELOPED, AND AN EXPERIMENT WAS RUN TO TEST 
THE MODEL. 

CORRELATION WAS SIGNIFICANT AT THE 99% LEVEL FOR SEVENTEEN RUNS ON 
THREE DIFFERENT GASES. HOWEVER A LEAST SQUARES REGRESSION ON THE DATA, 
PRODUCED A DESIRABLE SLOPE BUT A QUESTIONABLE INTERCEPT. AT WORST, THIS 
WOULD INDICATE THAT THE RANGEABILITY OF THE PREDICTOR IS GOOD ONLY AT 
HIGHER LEAKAGE RATES, BUT THE AUTHOR BELIEVES THAT THE DISCREPANCIES 
THAT OCCUR AT LOWER FLOWS ARE PROBABLY DUE TO ERRORS IN FLOWMETER 
CALIBRATION. 

IT WAS RECOMMENDED THAT A LARGER NUMBER OF TESTS OVER A WIDER 
VARIETY OF CONDITIONS BE RUN, AND THAT SOME SLIGHT MODIFICATIONS IN THE 
TESTING PROCEDURE BE MADE. THE AUTHOR BELIEVES THAT THE MODEL OR SOME 
MINOR VARIATION THEREOF WOULD BE AN ADEQUATE PREDICTOR FOR LEAK 
DETECTION. 


311 


IV 


TABLE OF CONTENTS 


I 

n 

BI 

IV 

V 

VI 


ACKNOWLEDGMENTS 

ABSTRACT 

SUMMARY 

TABLE OF CONTENTS 
NOMENCLATURE 
BODY OF TEXT 

6.1 INTRODUCTION 

6.2 MAIN TEXT 

6.2.1 DESCRIPTIVE INFORMATION 


6.2.3 RESULTS AND DISCUSSION 
vn CONCLUDING SECTION 

7.1 CONCLUDING REMARKS 
VIE APPENDIX 

LX REFERENCES 



312 



V 

NOMENCLATURE 


SYMBOL 

DEFINITION 

UNITS 

A 

ACOUSTIC VELOCITY 

M/S 

C 

SPECIFIC HEAT 

J/KG/K 

CD 

DISCHARGE COEFFICIENT 

LESS 

F 

SCALING FACTOR 

LESS 

G 

MASS VELOCITY 

KG/M**2/S 

H 

SPECIFIC ENTHALPY 

J/KG 

K 

SPECIFIC HEAT RATIO 

LESS 

M 

MOLECULAR MASS 

AMU 

MDOT 

MOLAR FLOW 

MOL/S 

P 

ABSOLUTE PRESSURE 

N/M**2 

R 

GAS CONSTANT 

J/MOL/K 

RHO 

DENSITY 

KG/M**3 

S 

CROSS-SECTIONAL AREA 

M**2 

T 

ABSOLUTE TEMPERATURE 

K 

V 

VELOCITY 

M/S 

SUBSCRIPT 

EXPLANATION 


a 

UPSTREAM CONDITION 


b 

DOWNSTREAM CONDITION 


c 

CRITICAL FLOW CONDITION 


gas 

GAS 


liq 

LIQUID 


o 

STAGNATION CONDITION 


p 

CONSTANT PRESSURE CONDITION 


* 

pro 

PROCESS CONDITION 


ref 

REFERENCE CONDITION 


s 

ISENTROPIC CONDITION 


V 

CONSTANT VOLUME CONDITION 
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VI 


BODY OF TEXT 


6.1 INTRODUCTION 

TOE LEAKING CONFIGURATION WAS ASSUMED TO BE THAT OF A SHARP EDGED 
ORIFICE. TOE CASES OF: 1) TOE SINGLE PHASE GAS, 2) TWO-PHASE FLOW, AND 3) 


ANY SINGLE PHASE FLUID (LIQUID OR GAS) MAY BE CONSIDERED AS THE PROCESS 
FLUID. THE ORIFICE EQUATION FOR A SHARP EDGEp ORIFICE HAYING A SMALL BETA 
IS APPLIED IN BOTH SITUATIONS AND TOE MOLAR FLOW OF THE REFERENCE FLUID 
IS DYNAMICALLY SCALED TO PREDICT THE MOLAR FLOW OF THE PROCESS FLUID. 

UNDER EACH CONDITION ONE HOPES THAT THE DISCHARGE COEFFICIENTS ARE THE 
SAME. IN REALITY THESE COEFFICIENTS HAVE A MAXIMUM RANGE OF ROUGHLY 
0.61 TO 1.0. THE LOWER EXTREME IS FAVORED FOR SHARP EDGED ORIFICES AT 
HIGH REYNOLDS NUMBER INCOMPRESSIBiE FLOW. TOE UPPER FOR CAREFULLY 
MACHINED VENTURI TUBES. INDUSTRIAL PRACTICE ROUTINELY ASSUMES A 
CONSTANT DISCHARGE COEFFICIENT IN TOE SIZING OF CONTROL VALVES, SAFETY 
VALVES, AND RUPTURE DISCS (1). 

6.2 MAIN TEXT 

6.2.1 DESCRIPTIVE INFORMATION 

A LABORATORY TEST WAS DESIGNED TO TEST FOR THE LEAKAGE OF VARIOUS 
GASES UNDER DIFFERENT CONDITIONS. THE GASES WERE FED THROUGH A PRES- 
SURE REGULATOR TO A BALLAST TANK. ATTACHED TO TOE BALLAST TANK WAS A 
SMALL BAR STOCK GATE VALVE THAT WAS SUGHTLY CRACKED. A SMALL 
ROTAMETER WAS THEN ATTACHED TO TOE OTHER END OF THAT VALVE. TOE RUNS 
CONSISTED OF VARYING THE UPSTREAM PRESSURE, FOR TOE VARIOUS GASES AND 
RECORDING TOE PRESSURE, TEMPERATURE, AND FLOW READINGS FOR EACH 
INDIVIDUAL RUN. SEVENTEEN DIFFERENT CONDITIONS WERE RECORDED. THREE 
DIFFERENT GASES WERE USED, NAMELY HELIUM, NTTROGEN, AND ARGON . 


THEN: 

dH + V*dV = 0, 

(1) 


P/RHO = 1000*R*T/M, 

(2) 


dH = Cp*dT, 

(3) 


Cp = Cv + 1000 * R/M 

(4) 


K = Cp/Cv, 

(5) 

& 

A**2 = [(partial of P)/(partial of RHO)]s 

(6) 

COMBINING THE ABOVE WE GET: 

P/RHO**K = Po/RHOo**K, (FOR ISENTROPIC STAGNATION) 

(7) 

& 

T/P**((K-1)/K) =To/Po**((K- 1 )/K) 

(8) 

SINCE 

G = RHO*V,& Gc = RHO*A 

(9) 

AND 

A**2 = 1000*K*R*T/M 

(10) 


Gc = RHO*sqrt(1000*K*R*T/M) 

(11) 

OR 

Gc = M *P/1 000/R/T*sqrt( 1 000*K*R*T/M) 

(12) 

Gc = P*sqrt(K*M/1000/R/T) 

ADJUSTING FOR STAGNATION 

Gc = Po* sqrt(K* (2/(K+ 1 ))* * ((K+ 1 )/(K- 1 ))*M/1 000/R/To) 

(13) 

(14) 

IF THE PROCESS IS NOT ISENTROPIC 

Gc = CD*Po*sqrt(K*(2/(K+l))**((K+l)/(K-l))*M/1000/R/To) 

(15) 


WHERE CD IS THE DISCHARGE COEFFICIENT. 

ACCORDING TO SHAPIRO (2)(CH.4,FIG.4.17)THE RANGE FOR CD IN CRITICAL FLOW 
VARIES BETWEEN 0.74 AND 0.85, AND DEPENDS ONLY ON THE RATIO OF DOWN- 
STREAM TO UPSTREAM PRESSURES. 

NOW LET US CONSIDER THE ALTERNATIVE CONDITION OF AN INCOMPRESSIBLE 
LIQUID UNDERGOING ISOTHERMAL FLOW IN A PERFECT CONVERGING-DIVERGING 
NOZZLE. THE FLOW IS AGAIN ISENTROPIC WITH STAGNATION UPSTREAM. THEN 
THE BERNOULLI EQUATION STATES: 

(Pb-Pa)/RHO + (Vb**2-Va**2)/2 = 0 (16) 

FOR Va = 0, LET Vb = V AND WE GET 

(Pb-Pa)/RHO + V**2/2 = 0 (17) 
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SINCE 

Gliq = RHO*V 

(18) 


(Pb-Pa)/RHO + (Gliq/RHO)**2/2 = 0 

09) 

OR 

Gliq = sqrt(2*(Pa-Pb)*RHO) 

(20) 

SIMILARLY FOR NON-ISENTROPIC FLOW 

Gliq = CD*sqrt(2*(Pa-Pb)*RHO) 

(21) 


WHERE CD IS THE DISCHARGE COEFFICIENT. 

DISCHARGE COEFFICIENTS FOR LIQUIDS FOR ORIFICES WITH SMALL BETAS RUN 
ABOUT 0.61 (3)(P.5-15).THERE IS NO RECOVERY OF THE LOSS (3)(P.5-17). 

THE THIRD POSSIBILITY WE WERE GOING TO CONSIDER WAS THE CASE OF TWO- 
PHASE CRITICAL FLOW. THE LITERATURE PROVIDES A LARGE AMOUNT OF INTER- 
ESTING APPROACHES TO THE PROBLEM, AND THE MODELING IS NOT THE OBVIOUS 
(4), (5), (6), (7), (8), (9). HOWEVER, THERE ARE INDICATIONS (4) THAT FOR SHARP 
EDGED ORIFICES WITH SATURATED LIQUID OR SATURATED VAPOR UPSTREAM THAT 
ONLY THE 100% LIQUID OR THE 100% VAPOR CASE NEED TO BE CONSIDERED. 
THEREFORE, FOR THE PURPOSES OF THIS REPORT ONLY SINGLE PHASE FLOW WAS 
INVESTIGATED. 

THEN FOR A GAS LEAK CONSIDER 

Gc = CD*Po*sqrt(K*(2/(K+l))**((K+l)/(K-l))*M/1000/R/To) (15) 

OR [Gc]pro [CD]pro*[Po*sqrt(K*(2/(K+l))**((K+l)/(K-l))*M/To)]pro 


(Gc]ref [CD] ref *[Po*sqrt(K*(2/(K+l))**((K+l)/(K-l))*M/To)]ref 

( 22 ) 


OR [Gc]pro [CD]pro 

* Fgas (23) 

[Gc]ref [CD]ref 

WHERE Fgas = [Po*sqrt(K*(2/(K+l))**((K+l)/(K-l))*M/To)]pro (24) 


[Po*sqrt(K*(2/(K+l))**((K+l)/(K-l))*M/To)]ref 

NOTING THAT Fgas IS A DIMENSIONLESS FACTOR D EPEND ING ONLY ON THE 
UPSTREAM PRESSURES, TEMPERATURES, AND IDENTITIES OF THE PROCESS AND 
REFERENCE GASES, OUR FOCUS SHIFTS TO THE RATIO OF THE DISCHARGE 
COEFFICIENTS - [CD]pro/[CD]ref. RE SHAPIRO (2) THE EXTREMES OF CD IN THE 
CRITICAL PRESSURE RATIO RANGE ARE FROM 0.74 TO 0.85 IMPLYING A MAXIMUM 
DIFFERENCE OF ABOUT 15% BETWEEN REFERENCE AND PROCESS CONDITIONS. 
HOWEVER IT IS EXPECTED THAT THE PRESSURE RATIOS WOULD BE SUFFICIENTLY 



SIMILAR OVER BOTH TEST AND REFERENCE CONDITIONS SO THAT THE RATIO 
[CDjpro/[CD] ref COULD BE TAKEN AT UNITY, AND 

[Gcjpro = [Gcjref * Fgas (25) 


Gliq = CD*sqrt(2*(Pa-Pa)*RHO) (26) 

[Gcjpro [CDJpre 

■— *. * Fliq (27) 

[Gcjref [CDJref 


WHERE Fliq = [sqrt(2*(Pa-Fb)*RHO)]pro (28) 


[Po*sqrt(K*(2/(K4T))**((K+l)/(K-l))*M/1000/R/To)]ref 



RELIABLY (3)(CH.5, FIG.5-20). 


HOWEVER, OTHER CONFIGURATIONS DEMONSTRATE CD’S BETWEEN 0.61 AND 1.0. 
CONSIDERING THE PREVIOUSLY MENTIONED GAS CD’S RANGING FROM 0.74 TO 0.85 
WE MAY ESTIMATE A RATE NO MORE THAN 20% LOWER THAN AN ESTIMATE BASED 
ON Fliq ALONE. NEVERTHELESS THIS IS AN EXTREME CONDITION, CONFIGURATIONS 
OTHER THAN SHARP EDGED ORIFICES TEND TO HAVE CD’S GREATER THAN 0.61, 
AND UNTIL EXPERIMENTAL DATA SHOWS OTHERWISE, ESTIMATES WILL BASED ON 
Fliq ONLY. THAT IS, THE [CDJpro/jCDjref RATIO WILL BE TAKEN AT UNITY. 


THEREFORE FOR GASES 

[Gcjpro = [Gcjref* Fgas (29) 

AND FOR LIQUIDS 

[Gliqjpro = [Gcjref * Fliq (30) 

SINCE MDOT = G*S/M (31) 

WE GET [MDOTJpro * [MDOTjref /(tMref]/[MproJ) * Fgas (32) 

OR [MDOTJpro = [MDOTjref /([MrefJ/[Mprol) * Fliq (33) 
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0.2.3 RESULTS AND DISCUSSION 


THE ROTAMETER DETERMINED FLOW RATES WERE CORRECTED FROM HELIUM TO 
THE PROCESS GAS INVOLVED (10)(FIG 8-27). THE DATA WERE THEN EVALUATED 
PER PROGRAM LEAX 0 DEVELOPED BY THE AUTHOR FOR THIS PROJECT. PICTORIAL 
REPRESENTATION OF THE RESULTS WERE THEN OBTAINED VIA* AN ADAPTATION OF 
THE IN HOUSE REGIS SOFTWARE. A LINEAR REGRESSION WAS APPLIED TO THE 
DATA AND A LEAST SQUARES FIT INDICATED A SLOPE OF ABOUT 0.9 AND AN 
INTERCEPT OF ABOUT 15 SCCM WITH SIGNIFICANCE AT THE 99% LEVEL. IDEALLY 
ONE WOULD EXPECT A SLOPE OF 1 .0 AND AN INTERCEPT OF ZERO. THE AUTHOR 
REGARPED THE SLOPE AS INDICATIVE OF THE RATIO OF THE CD’S WHICH WAS 
TAKEN AT UNITY. THOUGH THE 0.9 SLOPE IS QUITE SATISFACTORY, THE INTER- 
CEPT OF 15 SCCM WAS CONSIDERED TO BE HIGH AT LOW LEAK RATES. IN ORDER 



WTIH LEAX2 WHICH CORRECTED FOR THIS DISCREPANCY AND VERY SIMILAR 
RESULTS WERE OBTAINED. (TABLE 1, FIGURE 1, TABLE 2, FIGURE 2) 
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OBSERVED VALUE 


CONCLUDING SECTION 


7.1 CONCLUDING REMARKS 

X (PREDICTED) = 0.9 * X (OBSERVED) + B WAS SIGNIFICANT AT A LEVEL OF 99%, X 
REFERRING TO TOE MOLAR FLOW RATE. OF COURSE ONE WOULD IDEALLY EXPECT 
TOATX (PREDICTED) «X (OBSERVED). HOWEVER DUE TO TOE ROUGH NATURE OF 
”TE PREDICTOR 0.9 SEEMS TO BE ACCEPTABLE AS TOE SLOPE AND IS 
OF THE ASSUMPTION THAT TOE RATIO OF THE DISCHARGE COEFFICIENTS BE 
UNITY. TOE INTERCEPT B WOULD EXERT LITTLE INFLUENCE AT TOE HIGHER 
FLOWS BUT WOULD BE SIGNIFICANT AT LOWER FLOW RATES 


11113 izafcWfl m aas a nil. 
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THE RATIO OF TOE DISCHARGE COEFFICIENTS. VERY SIMILAR RESULTS WERE 
ACHIEVED, INDICATING THAT TOE RESIDUAL INTERCEPT MAY HAVE ARISEN FROM 
FLOWMETER CALIBRATION, IT IS RECOMMENDED, THAT TOE DATA RE RETAKEN ON 
A "TTME-WEIGH" BASIS TO DETERMINE WHETHER OR NOT TOE FLOW MEASURE- 
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APPENDIX 


LABORATORY DOCUMENTATION 
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